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We apply the formalism of the continuous time random walk to the study of financial data. 
The entire distribution of prices can be obtained once two auxiliary densities are known. These 
are the probability densities for the pausing time between successive jumps and the corresponding 
probability density for the magnitude of a jump. We have applied the formalism to data on the US 
dollar/Deutsche Mark future exchange, finding good agreement between theory and the observed 
CN ! data. 
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INTRODUCTION 



The continuous time random walk (CTRW), first introduced by Montroll and Weiss in 1965, Q, has a large number 
of applications to the modelling of many physical phenomena, particularly in the field of transport in disordered 
media J|, [| . In contrast to the standard random walk in which steps are made periodically, the CTRW is based on 
the assumption that the times between steps is random. The CTRW has been applied in many different fieds. These 
range from transport in amorphous materials 10], transport in turbid mediaj |, ^) , random networks , self-organized 
criticality ||, liquids |)|, electron tunneling jl0|], theoretical mechanics and earthquake modelling just to 
name a few. 

In this paper apply CTRW formalism to a phenomenon more related to social sciences than to natural sciences: the 
distribution of speculative prices. The first analytical approach to this class of problems was proposed and analyzed 
by Bachelier who in 1900 modelled stock price movements as an ordinary random walk where prices can go up and 
down, at fixed times, due to a variety of many independent random causes. This approach necessarily leads to the 



■ conclusion that the probability distribution of speculative prices is Gaussian 13 . In 1959 Osborne realized that, since 
stock prices are necessarily positive, it would be more convenient to consider returns instead of market values [ Ej . 
Thus, if S(t) is an speculative price (or the value of an index) at time t and 

2 : = HS(t)/s(o)\ (i) 

is the return up to time t. Then Z(t) is a Gaussian variable and S(t) is a log-normal process. Nevertheless, as Kendall 
first noticed in 1953 the normal density fits financial data very poorly in the tails of the distribution. As an 
example: the probability of an event corresponding to five or more standard deviations is up to 10 4 times larger than 
the one predicted by the Gaussian distribution. Therefore, empirical price distributions are highly leptokurtic. The 
existence of these "fat tails" was precisely what led to Mandelbrot in 1963 to propose the Levy distribution for stock 
market prices (l(|. There is, however, a drawback to this approach: no finite moments exist beyond the first and 
this is certainly a severe limitation of the model. Moreover, the Levy distribution has been tested against data in 
. many situations, always with the same conclusion: the tails are far too long compared with actual data. In any case, 
as Mantegna and Stanley have recently shown p7| , the Levy distribution fits very well to the center of empirical 
distributions — surprisingly much better than the Gaussian density — and it also shares the scaling behavior that 
appears in data. 

Recently a new market model was proposed to fill the gap between Gaussian and Levy distributions Jlc| . The 
model, which was based on a continuous superposition of jump processes, explains the appearance of fat tails and 
self-scaling but still keeps all moments finite. It reproduces price distributions quite exactly, particularly those of 
tic-by-tic data. In this paper we want to address the problem from a different point of view. Thus we assume that the 
evolution of prices can be modeled by a CTRW. This allows us to calculate the distribution of speculative prices. The 
paper is organized as follows. In Sections II and III we set the general formalism and derive the exact distribution of 
prices and volatility. In Section IV we derive some asymptotic results mostly valid for long times. In Section V we 
apply the model to real data, namely the US dollar/Deutsche Mark future market. Conclusions are drawn in Section 
VI. 
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GENERAL FORMALISM. THE DISTRIBUTION OF PRICES 



Wc define the zero- mean return X(t) by 

X(t)=Z(t)-{Z(t)), (2) 

where Z{t) is given by Eq. (|l|) and (Z(t)) is its average. We now suppose that X(t) can be described in terms of a 
CTRW. In this picture X(t) changes at random times to, t%, t 2 , ■ ■ ■ , t n , • ■ ■ and we assume that the intervals between 
successive steps (which we call "sojourns") T n = t n — t n —i (n = 1, 2, 3, • • •) are independent and identically distributed 
random variables with a probability density function given by tp(t), i.e., 

ip(t)dt = Prob{i < T n < t + dt}. 

At a given sojourn the zero- mean return X{t) undergoes a random change giving rise to the random variable AX n = 
X(t n ) — X(t n —i) which is described by a probability density function defined by 

h(x)dx = Probjx < AX n < x + dx}. 

In this formulation of the problem we choose a function p(x, t) to be the fundamental function, where p(x, t)dxdt is 
the joint probability that an increment in return, X(t), is added whose magnitude is between x and x + dx and that 
the time between successive turns is between t and t + dt. The condition that there is no net drift will be assured by 
requiring that p(x.t) is an even function of x. We can form two marginal densities out of p{x,i): the pausing-time 
density i/j(t) for the time between successive pulses 



$(t) = / p(x,t)dx, (3) 

J — OO 

and the probability density function (pdf) for the changes in a single jump, h(x), where 

h{x) = / p(x,t)dt. (4) 



The value of the return X(t) at time t will by given by the random value of the height at t (see Fig. [j]). We are 
interested in the probability density function of this variable, p(x,t). 

We also assume that between successive steps the time evolution of X(t) is linear (see Fig. |l|). The assumption 
of linearity between jumps is arbitrary. We could have used step functions instead, in such a case the return would 
evolve discontinuously and during any sojourn the value of the return (and hence the price) is that of the last jump. 
On the other hand, the linear choice for the return X(t) implies an exponential growth in price S(t) (cf Eq. (pi)). 
This exponential behavior is an inherent feature of any financial market. 

Apparently the choice of linearity is against causality because one has to know future return values, X(t n ), in 
order to draw the segment joining X(t n -i) and X(t n ). Fortunately for the soundness of our model this is not true. 
Indeed, notice that our process is completely equivalent to a random process evolving linearly with a slope changing 
at the random times t n (n — 1, 2, 3, • • •). If the slope of the segment joining X(t n _x) with X(t n ) is independent on 
tmin+ii • ' ■) then the random process is causal. Moreover, if the slope is independent on its previous values then the 
process is also Markovian. 

Let us calculate the form of p(x, t) prior to the first jump. This function will be denoted by po(x, t) and, due to the 
linear evolution of X(t) between steps, it reads 



M) 



p (x, t)= dT p(y, T)8[x- y -) dy. (5) 



where we have assumed that the initial jump occurred at t = 0. In terms of po and p we have that the pdf p(x, t) for 
the return at time t is given by 

ft />oo 

p(x,t) = po(x,t) + / dt' p(x',t')p(x — x',t — t'). (6) 

JO J -oo 

This equation has been derived from the consideration that at time t, the process is either within the very first sojourn, 
this given by the first term on the rhs of Eq. (^|), or else the first sojourn ended at time t' < t, at which time the 
return had value x' , and from [x 1 t') the process was renewed. 
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It is possible to solve Eq. (Q) by means of a joint Fourier-Laplace transform. To this end let us denote by 

/>oo />oo 

p(uj,s)= dte- st du;e iwx p(x,t), 







the joint Fourier-Laplace transform of p(x, t). Then the convolution theorems applied to Eq. yield 

po(u),s) 

P(u,s) = - -, (7) 

1 - p(w,s) 

where Po(lu, s) and p(uj, s) are respectively the joint Fourier-Laplace transforms of functions paix, t) and p(x, t). Recall 
that po{x,t) only depends on p{x,t), thus the return pdf p{x,t) is exclusively determined by the form of p(x,t). 

Unfortunately the form of p(x, t) is very difficult to determine from the available data. More easily accessible are 
the marginal densities of p(x, t) in Eqs. (||) and (Q). It is therefore essential to assume a functional relation between 
p(x,t) and its marginal densities ip{t) and h(x). The simplest choice would be based on the assumption that return 
increments and their duration time are independent random variables. In this case 

p(x,t) = h(x)i/>(t). (8) 

However, this assumption does not seem to be realistic since one certainly expects some degree of correlation between 
return increments and their duration, while Eq. (^) implies complete independence between increments and sojourn 
times. Following that intuition, for the rest of the paper we will mostly assume that the density p(x, t) is such that 
its characteristic function p(w,t) has the functional form: 



p(u,t) = ip 



t 



h(oj) 



(9) 



where ip(t) is the pausing time density, which we assume to be a decreasing function of time, |19 , and h(uo) is the 
characteristic function of h{x). Notice that Eq. (S) has been chosen to satisfy Eqs. (^|) and (Q). 

As we have mentioned, assumption given by Eq. (^|) is at least intuitively plausible because it implies that one 
must wait for a long time in order for a large variation of return to occur. In other words, major increments of the 
return are very infrequent. We will prove this by showing that sojourn time T = t n — t n _i and return quadratic 
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increments AX 2 = \X(t n ) — X(t n _i)] 2 have a positive correlation, then increasing return variations imply increasing 
sojourn times and vice versa. In effect, we define the correlation function between AX 2 and T by 

r = (AX 2 T) - (AX 2 )(AT). 

We can easily evaluate the cross average (AX 2 T) using the joint characteristic function /5(w, i).Thus 

duo 2 



(^ 2 T) = - 



tp(oj, t)dt 



;=0 



which after using Eq. (Q) yields 



Hence 



(AX 2 T) = 2(AX 2 )(AT). 
r = {AX 2 ){AT) > 0, 

as we meant to prove. 

Observe that assumption (|^) allows us to write the joint Fourier-Laplace transform of pix. t) in the form 

p((jj, s) = sh(uj) 



(10) 



where "0(s) is the Laplace transform of the pausing-time density ip(t). Likewise, from Eq. (||) we see that the 
transformed density fio(uj, s) can be written as 

d f 1 dz 



6»a 



-p(ujz,sz), 



which, after using Eq. @, reads 



Po(uj,s) = — / h 2 {ojz)ip' szh{u}z) 



dz, 



(11) 



(12) 



where the prime denotes a derivative. In terms of the transformed densities ip(s) and h(oj) the formal solution to the 
problem given by Eq. (0) can be written in the following more explicit form 



p(w,s) = 





szhiuiz) 


dz 


1 - h(u)ip 


sh(uj) 





(13) 



Subject to the assumption in Eq. Eq. (13) furnishes a complete solution to the problem and it can be a convenient 
starting point for numerical methods when further analytical insight is unavailable. 

We finish this section with an example. Suppose that the random times to,ti,t 2 , ■ ■ ■ in which the return suffers 
random increments form a Poisson set of events, suppose also that these random increments are distributed according 
to a Laplace density. Then densities ip(t) an d h(x) are respectively given by 



ip(t) = Ac 



-XI 



(14) 



where A 1 = (T) is the mean sojourn time and 7 > is such that (AX 2 ) = 2/^ 2 is the jump variance. In this case 
h{uS) = 1/(1 +u! 2 /j 2 ) and the characteristic function of the joint density p(x,t) is given by 

p( w ,i)=Aexp{-A(l+^ 2 /7 2 )i}, (15) 

which is the convolution of a Poissonian density Xe~ xt and a Gaussian density with zero mean and variance 2Xt/j 2 . 
Finally, the inverse Laplace transform of the formal solution in Eq. ( |l3"| ) now results in an explicit expression for the 
characteristic function of the problem given by the time convolution: 



p(cj,t) = po(u>,t) + A 



po{u>,t - t')dt', 



where 



Po(u,t) = A 



exp 



-At 



u z t 2 

J 2 T 



dr. 



(16) 



(17) 



Note that numerical analyzing Eqs. (|I1|)-(|l7j) is straightforward. 
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THE VOLATILITY 



Besides the pdf p(x,t), which provides all possible information about the problem, there is another quantity of 
considerable practical interest: the return variance. In our analysis this quantity, called "volatility" in the terminology 
of finance, has the advantage that it does not require the knowledge of the entire jump distribution h(x). It suffices 
to know the pdf tp(t) and the first two moments of h{x). 

Let us denote by (X n (s)) the Laplace transform of the nth moment of the return: 

(X n {s)) = / e- st (X n (t))dt. 



(18) 



This can be written in terms of the joint Fourier-Laplace transform of p(x, t) by 

-n d n p(uj,s) 



duj n 



Since we assume that there is no net drift in the evolution of X(t), i.e., p(—x,t) = p(x,t), this means that all odd 
moments associated with p(x, t) are zero. That is, 

(i? 2 "- 1 ^)} -0, (n = 1,2,3,.-.), (19) 

where 

/oo 

x m p(x,t)dx (20) 
-oo 

is the m-th moment of the return increment. All of this means that random jumps during any sojourn are unbiased 

and, in particular, that their average is zero. Note that since the jump density h(x) is a marginal density of p(x,t), 
condition ([l9]) implies that 

hf- 2n ~ 1 \0)=0, and pi n = {~l) n U 2n \0), (21) 

(n = 1, 2, 3, • ■ •), where h^ m \0) is the m-th derivative of the characteristic function, h{ui), of the jump density. Notice 
that another direct consequence of the unbiased assumption given by Eq. ( |l9| ) is that all odd moments of the return 
process vanish: 

(A 2 »- 1 (i))=0, (n = 1,2,3,..-). (22) 
Starting from Eq. (|l|) and using Eq. (0) and Eqs. we obtain 

<**(.)> = <5 2 (')> + <*(«)>, (23) 

i - i*{s) 



where (R 2 (s)} and ip(t) are the Laplace transforms of (R 2 (t)) and ip(t) respectively, and 

( x o( s )) = - 



The substitution of Eq. ([ll]) into Eq. (Q) and some simple manipulations finally yield 

(Xl(s)) = --(R 2 (s)) + - f z(3( S z))dz, 

S S In 



(24) 



(X 2 (s)) = 2/ f / z(R 2 (sz))dz. (25) 



and Eq. (|23j) implies 

(X 2 (s)) = 

1 - tp(s) Jo 

As we have mentioned, the independent model given by Eq. (JsJ) cannot be used for describing actual markets, but 
for the sake of completeness we will also give the general expression of the volatility associated with the model. This 
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result will serve to illustrate an important point regarding the asymptotic behavior of the volatility which will be 
discussed nearly at the end of the next section. 
Using Eq. (pi) we have 



and 



(R?(s)) =/j 2 V(s), 



(X 2 (s)) = 2/i2 / S [ zip(sz)dz, 



1 - ip(s) Jo 

while for the dependent model exemplified by Eq. and after using Eq. (|l0|) we have 



(26) 



(R 2 (s))= m - 8f/>{a) 



Then 



zip(sz)dz 



(27) 



Observe that both Eqs. ( p6| ) and (]27|) depend only on the pausing-time density and the second moment of the jump 
/i2- In this case and for the example given in Eqs. ([T^)-(|l7|) we can explicitly write 



(x 2 (t)) 



2Xt 



3^2 ry2 3^,2 



(-1 + 2Xt + X 2 t 2 )e- Xt 

+ \\ 2 t 2 (l + \t/3)E 1 {Xt), 



(28) 



where Ex (At) is the exponential integral. Using both the short time behavior and the asymptotic behavior of Ei(Ai) 
pot we can easily see that 



and 



(X 2 (t)) ~ -(2/ 7 2 )A 2 i 2 lnAi, (t < A" 1 ), 



(X 2 (t)) ~ (2A/ 7 2 )t, (i»A- 1 ). 



In this particular case we thus observe an anomalous diffusion-like behavior at short times and a diffusion-like behavior 
at long times. We will see next that this is a general feature of the model rather than a peculiarity of the example 
given by Eq. ( O ) . 

ASYMPTOTIC RESULTS 

In this section we obtain some approximate results regarding the characteristic function and the volatility of the 
process. These results will mostly refer to the behavior of the probability distribution for large t and large x. To this 
end we will assume that the Laplace transform of the pausing-time density, tp(s), has the following series expansion 



b n s 



a-\-n— 1 



(29) 



where a > 1 is a non integer number. Note that Eq. (|29|) is a fairly general assumption, because when b n = for all 
n then all moments (T n ) of ip(t) exist. In such a case a n = (—l) n (T n )/n\. On the other hand if ip(t) only possesses 
N > 1 moments, then N +1 < a < N + 2 and a n = (-1)" (T n ) jn\ only for n < N. Using Eq. Q we see that p(w, s) 
given by Eq. @ reads 



p(w,s) = h{u) 



71=1 



n=0 



(30) 
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while po(uj, s) is (c/Eq. (|l 



Pi) 



{w,s) = -J2a n s n V«M - J2(a + Ti)b n s a+n 2 </> Q +n-i (v), 



(31) 



where 



</>/c(w) 



z* _1 ft k (w«)<fo. 



(32) 



We want to obtain an asymptotic expansion of the pdf p(x, t) valid for large t and |x|. As is well known, the large 
t behavior is equivalent to the small s behavior of the Laplace transform. Similarly the large \x\ behavior correspond 
to the small u> behavior in the Fourier domain. Having this in mind from Eq. (B(J) we have for a > 3: 



l-p(u,s) l-h{Lu)~ ai h 2 (uj)s + 0{s 2 ) 



or equivalently, 



1 - p(u, s) ai h 2 (w)s{l - [1 - h(uj)}/[ ai h 2 (u;)s] + 0(a)} 



(33) 



Recall that u> is also small and in this case h(u>) ~ 1. We now assume that, in spite of s being small, the range of 
small values of lo that we will consider is such that 



< 1. 



1 - h(w)\ 
\aih 2 (io)s\ 

We can thus expand the rhs of Eq. (^) with the result 

- — — [l + l -J^- + 0(s) 
l-p(w,s) aih 2 (w)s Y a x h 2 (u))s 

On the other hand from Eq. (|3l]) we can write 

Po(^,s) = -ai0i(w) + O(s), 

where <f>\ (lu) is given by Eq. @. The substitution of Eqs. @ and (||) into Eq. @) yields 



h 2 {uo)s 



1 - h(u) 
1 + - v ' + O(s) 
aih 2 (uj)s 



up to the leading order, Eq. (|36|) yields 



where a\ — ~{T) is equal to the negative of the mean time between successive jumps and 



(34) 



(35) 



(36) 



(37) 



(38) 



Thus by virtue of Tauberian theorems the asymptotic expression for large t of the characteristic function p{u,t) will 
be given by the Laplace inversion of Eq. (37): 



P( U ,t) ~ -N(w)-=r, (t -►«>). 



(39) 



In this asymptotic case the market volatility could be evaluated starting from Eq. ( p7| ) and then following the 
procedure just described which involves the use of Eq. (p9|). Nevertheless, it turns out to be much simpler to directly 
evaluate {X 2 (t)) using the asymptotic expression for p{w,t) given by Eq. (p59). We thus have 



(X 2 (t)) ~ m t/(T) (t oo). 



(40) 
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Therefore, in this regime the volatility grows linearly with time, which suggests a diffusion-like behavior of the model 
at long times. 

It is also interesting to calculate the behavior of volatility at short times. To this end we first note that as t — > no 
shift in return has occurred with high probability. Consequently the probability density function p(x, t) of the process 
is approximately given by po(x, t): 



p{x,t) ~p (x,t), 
and the approximate expression of the volatility is 



(t -> 0), 



which, after using Eq. (0) yields 



(X 2 (t)) 



xpo(x,t)dx, (t — > 0), 



where (R 2 (t)), the second moment of p(x, t), is defined by Eq. (|2C 
have 



(t -> 0), (41) 
In terms of the characteristic function p(uj,t) we 



(R 2 (t)) = 



d 2 p(L0,t) 



and from Eq. (|9|) we get (see Eq. 



(R 2 (t)) = -p 2 t^'{t). 



Substituting this equation into Eq. (41) yields 

(X 2 (t)) ~ -p 2 t 2 



new) 



(i-0), 



lt/(T) £ 

where we have used dimensionless units in writing the integral. An integration by parts yields 

r tWA d ^^ [t)Ht/{T)) 

Jt/(T) S Jt/(T) 



where we have assumed that i/j'(t) decreases fast enough at infinity. It is easy to convince oneself that as t 
dominant term on the right hand side of this equation is the first one. Hence, 



the 



t/(T) 



d£~-^(t)M*/CO)- 



Finally, 



(X 2 (t))^^'(0)< 2 ln(V(T)), (t<<T». 



(42) 



Since t 2 \ ln(t/(T))\ < t when t <C (T), the volatility grows slower than normal diffusion at short times. Therefore, 
the model exhibits an anomalous diffusion- like behavior at short times. This peculiar behavior of the volatility, i.e., 
anomalous diffusion at short times and ordinary diffusion at long times, is a characteristic feature of the model, and 
we will see next that there seems to be empirical evidence of such a behavior in real markets. 

It could be argued that the anomalous behavior of the volatility at short times is a spurious consequence of the 
form of po(x,t) which, in turn, is the result of the arbitrary assumption that the time evolution of (X 2 (t)) between 
successive steps is linear (cf Eq. (||) and Fig. [I]). We will prove that this is not the case and that the anomalous 
behavior of volatility is, at least, a direct consequence of assumption (^J) which relates return increments with their 
duration. In effect, suppose that the expression for po(x,t) given by Eq. (|^) is valid but that return increments and 
time intervals are independent random variables. Then from Eq. (H) we have 



p(w,t) = h(uj)ip(t). 



(43) 
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Therefore (R 2 (t)) = ^{t) and Eq. @ reads 

(X 2 (t))c^ 2 t r^lldr. (f-0), 

For t sufficiently small we may write 

(X 2 (t)) ~ /i 2 ^(0)t, (i^O), 

and the model presents a diffusion-like behavior at short time scales. Hence, it is the dependence between return 
increments and their duration the reason for the anomalous behavior of the model at short times. We finally note 
that this provides a test for the validity of assumption in Eq. (||) since if actual data do not support diffusion-like 
behavior at short times then the assumption of independence between increments and their duration is inaccurate. 

We finish this section obtaining the asymptotic long time behavior of the volatility for the independent model given 
by Eq. (||). In this case (X 2 (s)) is given by Eq. ( p6|) which, for small s and after using the expansion (p9|), reads 

\-l(T)s + 0( S 2 ) . 

Whence 

(X 2 (t)) ~» 2 t/(T), (t-oo) 

and, for the independent model, the volatility also has the same diffusion-like behavior at long time as that of the 
dependent model (c/Eq. (ff"o|)). 



SPECIFIC RESULTS 



We will now apply the formalism presented in previous sections to analyze the distribution of returns corresponding 
to the future US dollar/Deutsche Mark exchange (tic-by-tic data from 1993 to 1997) |2l], 

Before proceeding further we need to comment on a key point regarding the nature of markets. In a majority of 
works on the subject it is implicitly assumed that statistical properties of the economy are stationary all over the 
time. This is certainly inaccurate, for one does not expect current behavior of the market to be that of the market 
in, say, 1930. Here we presuppose a less restrictive assumption than that of complete stationarity: we suppose that 
markets are stationary over shorter periods of time, say one or two decades. Since we work on high frequency data 
and these data are only available since the early 1990's, our assumption of "local (in time) stationarity" seems to be 
consistent with the data. 

As explained above, our first task is to infer from these data which forms for ip(t) and h(x) are plausible. In Fig. 
U we plot the experimental pausing-time density i/)(t). We can see there that an excellent fit to the data is provided 
by the following pdf 

m = n+ At)« ' ( 3<a<4 )' ( 44 ) 

where a = 3.47 and A = 2.73 x 10~ 2 s -1 . Since 3 < a < 4 the first two moments of tp(t) are finite while the rest of 
moments do not exist. The Laplace transform of tp(t) is therefore of the class given by Eq. (p9f). The mean sojourn 
time and the second moment are 

(T) = — r (T 2 ) = 2 *"' (45) 
a — I (a — 2) (a — 6) 

For the dollar/mark future market the experimental mean sojourn time evaluated from data is (T) exp = 23.65 s, in 
satisfactory agreement with the theoretical prediction of (T) = 24.85 s evaluated from Eq. (|45|) and based on ansatz 

In Fig. ^ we plot the experimental jump density h(x). We see there that the experimental h(x) can be considered 
a symmetric function of the return increments x. A good fit is also given by a power law with a greater exponent 
than that of ip(t): 

h W = W + \*\hT (5</?<6) ' (46) 
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1 10 100 1000 

Time between sojourns (in seconds) 

FIG. 2: Empirical distribution of the time between transactions, corresponding to the operative of the closest-to-maturity 
Deutsche Mark future (in the US market). The analyzed data range from January, 1993 to December, 1997. The pdf of the 
sojourn times, if>(t), clearly follows a power law. The solid curve represents the fit we propose in the main text. 



where (3 = 5.52 and 7 = 2.64 x 10~ 4 . Again, the pdf in ( firf ) has its four first moments finite and the rest are infinite. 

Power law densities like (Q) and (ji^) have been recently suggested for describing several market models such as 
individual companies |p3| or market indices [2^] . 

We now have all the ingredients to obtain a complete analysis of the dollar/mark exchange market. Unfortunately 
the densities given by Eqs. ([5]) and ( f46| ) make very problematic the exact evaluation of p(x, t) by means of the 
Fourier-Laplace inversion of Eq. (113). We will use instead the approximate results obtained in Sect. IV, and in 



particular the asymptotic expression Eq. (39) for the characteristic function p(u,t) valid for large t. To this end we 



must have an expression for the jump characteristic function, 

/>oo 

h(uj) — 2 / h(x) cos uixdx. (47) 



Substituting Eq. ( flq ) into this equation enables us to obtain the exact h(u>) in terms of a combination of incomplete 
gamma functions of imaginary arguments. However that expression is clumsy for practical purposes and, since we 
are mainly interested in the behavior of the tails of p(x, t), we will use a simpler expression for h(u>) valid when u) is 
small. 

We therefore define 

g(w) = 1 - h(oj) 



then, using Eqs. (46)-(|47j) and taking into account both normalization and symmetry of h(x), we write 

/>oo 

g(oj) = 2 / [1 — h(x)] cos ojxdx. 



Substituting Eq. (j4^) into this equation and recalling that 5 < (3 < 6, after successive integrations by parts we find 

r(/3-3). , 2 r(/3-5). , 4 



sin x 

ax. 



r(/3-i) Jo xP-^il + oj^/xY 
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X 




Return variations between sojourns 

FIG. 3: Empirical distribution of the logarithmic changes between transactions, in the Deutsche Mark/US dollar future mar- 
ket. Positive variations (increments) and negative variations (decrements) exhibit approximately the same behavior, thereby 
supporting our assumption of the symmetry of h(x). The plot also suggests the presence of a power law, and it includes a 
graph showing the shape of h(x) in Eq. (|46|), using the parameters reported there. 



As to — > we make the approximation 



sini f°° sin a; 

rdx ~ / — = — -ax 



Hence 



o xP- 5 (l + uj 1 / x y- 5 J x?- 5 2r(/3-5)sin7r(/3-5)/2' 



r (/ ? - 3 ),...^2 , r(/3-5). . 4 



{^f-\ (48) 



2r(/3-l)sin7r(/3-5)/2 



The Fourier inversion of this approximation will give us the behavior of h(x) as x — > ±oo. Then, neglecting delta 
function terms (which obviously do not contribute for large values of \x\), we have 



h(x) ~ — ; — — — 7^ 1 / lu 13 1 cos Loxdco, 

{ ' 2r(/3-l)sin7r(/3-5)/2 1 J 

(x — > ±cx)).The integral appearing on the right hand side of this equation is convergent in the sense of generalized 
functions and reads j2(]] 

a i , r(/3)sin7r(/3-5)/2 
lo^- 1 cos uxduj = ^ — ^ — . (49) 

jo \ x r 

Therefore, 

and the tails of the jump distribution follow a power law with exponent 3. This is consistent with Eq. (^) which in 
turn proves the soundness of the approximation in Eq. (48). 
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FIG. 4: Empirical probability density function p(x,t) for a set of time lags t, ranging from 15 seconds to 2 minutes. The 
model leads to a power-law decay, governed by an exponent f3, in all cases. This exponent is precisely the one appearing in the 
power-law for the jump density h(x). Since we have chosen exponentially growing values for t, different times correspond to 
equally spaced tails and this indicates a linear growth with time. We only show tails for positive increments, tails for negative 
increments behave in the same way. 



Let us now prove that the tails of the entire distribution p(x, t) also obey a power law with the same exponent 
for any time t sufficiently large. Indeed, the substitution of Eq. ( fi"8| ) into Eq. (|3|) yields 

N{lj) ~ A/( 7 w) /3 -\ (51) 

where 

7T 

M = 2T{0- l)shi7r(/?-5)/2' (52) 

In writing Eq. ( |5l"| ) we have taken into account the fact mentioned above that integer powers of to do not affect the 
behavior at the tails. In this situation the asymptotic expression of p(u>,t) given by Eq. ( |39| ) can be written as 

p(w,*)~-Mt( 7 w)"- 1 /(T), (53) 
and the Fourier inversion of Eq. (|5^) finally reads ^6) 

(0 - DiV 3 " 1 

2(T) W ' (N ^°° ) - (M) 

Hence tails of p(x, t) decay following the same power law as that of the return increment distribution h(x). This 
prediction of the theoretical model is confirmed by actual data. In Fig. ^ we show the empirical p(x, t) for the 
dollar/mark future exchange and for different values of time t, ranging from 15 seconds to 2 minutes. The empirical 
distribution clearly shows, for all these times, a power-law decay with exponent w 5.5, which coincides with the 
decaying exponent of h(x) thus confirming the predictions of the CTRW model. Moreover, Eq. (|54|) predicts the 
linear growth of tails with time. This linear growth is indeed observed in Fig. ^ where different times correspond to 
properly spaced curves. 

Let us finally and briefly comment on the volatility. In Fig. || we plot the volatility (X 2 {t)) for the dollar/mark 
future market. We can see there that the experimental volatility clearly shows two different regimes: at short times we 
observe a subdiffusion-like behavior while at long times the volatility clearly appears to have a diffusion- like behavior. 
Both regimes are consistent with the CTRW dependent model. We also note from Fig. || that the transition between 
these two regimes occurs at t ~ 25 s which is approximately equal to the mean sojourn time (T). Again this transition 
from subdiffusion to diffusion around time t ~ (T) has been predicted by the dependent model. 
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10° 10 1 <T> 10 2 10 3 10 4 

Time (seconds) 

FIG. 5: The experimental volatility (dots) shows two different regimes. At short times there is a subdiffusion-like behavior (dot- 
ted line) while at long time the volatility grows linearly with time (solid line) . Transition between regimes occurs approximately 
at time t ~ (T). These facts are in agreement with theoretical predictions. 



CONCLUSIONS 

In this paper we have applied formalism based on the CTRW to the random movement of market prices. The 
formalism depends on the conjecture from data of two densities: the pausing-time density ip(t) and the jump density 
h(x). The assumption that both densities are independent necessarily leads to the conclusion that the volatility of 
the return process has a diffusion- like behavior, i.e. grows linearly with time, at any time scale. The hypothesis 
of independence is perhaps the simplest assumption one can make. However, it does not seem to be realistic since 
return variations and their duration are certainly correlated, at least in many markets. We have therefore proposed a 
dependent model in which large return increments are infrequent. With this assumption the model predicts that the 
volatility should behave in an anomalous diffusive way at short times, something that is seen in some markets. 

The CTRW formalism allows us to obtain a closed expression for the joint Fourier-Laplace of the entire process 
which constitutes a convenient starting point for numerical analysis when no further analytical manipulations can be 
made. We have also obtain an asymptotic long-time expression for the characteristic function of the return process, 
this expression is valid as long as the pausing-time density possesses a hnite first moment (T). 

We have applied the formalism in a study of the US dollar/Deutsche Mark future exchange market. We have used 
tic-by-tic data from 1993 to 1997. Data show that ip(t) and h(x) are very well described by power-law densities (c/ 
Eqs. ( fl4| ) and (f46|)). We have showed that (i) the tails of the return distribution follow a power-law with the same 
exponent as that of h(x). (ii) The volatility has a diffusion-like behavior at long times and an anomalous diffusion- like 
behavior at short times. Both conclusion agree with experimental data. 

Let us finally mention that we have not been able to apply the CTRW formalism to market indices, such as the 
S&P 500, since any index is an average of many prices and thus indices are recorded at fixed times. This contradicts 
the underlying assumption of the CTRW, i.e., that the time between successive changes is random. Therefore, the 
formalism presented herein is valid and applicable to single companies, currency exchange and commodities, while for 
market indices other formalisms, like the one presented in fljg|j , are necessary. 

This work has been supported in part by Direction General de Investigation under contract No. BFM2000-0795 
and by Generalitat de Catalunya under contract No. 2000 SGR-00023. We thank Marian Boguna for many comments 
and suggestions to improve the manuscript. 
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